Growth Inside a Corner: The Limiting Interface Shape 
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We investigate the growth of a crystal that is built by depositing cubes inside a corner. The 
interface of this crystal approaches a deterministic growing limiting shape in the long-time limit. 
Building on known results for the corresponding two-dimensional system and accounting for basic 
three-dimensional symmetries, we conjecture a governing equation for the evolution of the interface 
profile. We solve this equation analytically and find excellent agreement with simulations of the 
growth process. We also present a generalization to arbitrary spatial dimension. 
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Growing interfaces constitute a venerable subject, but 
the proper continuum framework to account for this 
growth was developed not so long ago [1] . A detailed and 
beautiful description of fluctuations of one- dimensional 
growing interfaces has been proposed [2j [3], culminat- 
ing in a recent solution of the KPZ equation [J. For 
real applications, two-dimensional growing interfaces are 
much more important, but their governing stochastic 
continuum equations [1 remain unsolved. Nevertheless, 
the analysis of two-dimensional growing interfaces is not 
hopeless. Indeed, although interface fluctuations have at- 
tracted the most attention, they become less important 
as the interface grows. The limiting shape — the aver- 
age interface profile in the long-time limit — is the more 
primal characteristic. 

If growth begins from a flat substrate, the interface ad- 
vances at a constant average speed, so only fluctuations 
matter. In numerous applications, however, the limit- 
ing shapes are curved and are known only in rare cases. 
One such example is the 2+1 dimensional Gates- Westcott 
model for vicinal interfaces, which was solved by a free- 
fermion mapping [5]. This growth process exhibits loga- 
rithmic height correlations and therefore does not belong 
to the strong-coupling KPZ universality class. Average 
interface profiles are also known for certain anisotropic 
2+1 dimensional growth models [6l[7]. However, even for 
the most basic isotropic growth models limiting shapes 
are not known. For example, for the two-dimensional 
Eden-Richardson model ^ the limiting shape is un- 
known, although the statistics of its fluctuations are un- 
derstood (and belong to the KPZ universality class). 

Here we investigate the limiting shape of a crystal that 
grows inside a corner. This process can be defined in ar- 
bitrary dimension and on any lattice (with an appropri- 
ately defined 'corner'). We specifically consider a cubic 
lattice, where the corner is the initially empty positive 
octant. Starting at t = 0, elemental cubes are deposited 
at unit rate onto inner corners (Fig. [T]). Initially, there 
is one inner corner and thus one place where a cube can 
deposit. After this first event, there are three available 
inner corners that can accommodate the next cube. The 



interface shape becomes smoother as it grows and ulti- 
mately approaches a deterministic limiting shape. 




FIG. 1: (color online) Upper left: 3d crystal of volume 4. 
The next elemental cube can be deposited at one of 6 inner 
corners. Right: Crystal at t = 140. 

The corner growth model admits a dual interpretation 
as the melting of a three-dimensional cubic crystal by 
erosion from the corner. There is also a magnetic inter- 
pretation in which plus spins are initially assigned to each 
site inside the corner and minus spins to exterior sites, 
with the spins endowed with zero-temperature Glauber 
spin- flip dynamics [9 in a weak negative magnetic field. 
This dynamics allows only plus spins at inner corners to 
flip and thus is isomorphic to the corner melting problem. 
The magnetic interpretation naturally suggests consider- 
ing the system in zero magnetic field, which results in a 
growing interface whose characteristic scale grows diffu- 
sively rather than ballistically. Other modifications in- 
volve changing the initial condition; e.g., depositing the 
cubes onto a planar substrate (the 'hypercube stacking 
model' [10 ) leads to a trivial limiting shape but is better- 
suited to studying non-trivial height fluctuations. 

In what follows, we use the language of deposition; 
most importantly, we allow only deposition events and 
no evaporation. Growth inside a two-dimensional cor- 
ner is well understood by mapping the corner growth 
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process onto the one- dimensional asymmetric exclusion 
process [11] ; fluctuations in this limiting shape have also 
been computed [121 [13]. three dimensions, the corner 
growth model can be mapped into an infinite set of cou- 
pled exclusion processes in the plane, also known as the 
'zigzag model' [l4l[T5! . Unfortunately, no exact solutions 
are known for such planar interacting particle processes. 

Here we focus on the limiting shape in three (and 
higher) dimensional corners. Our analysis relies heav- 
ily on insights gleaned from the limiting two-dimensional 
corner interface shape [11 . In two dimensions this lim- 
iting shape y{x;t) evolves according to the equation of 
motion [iQ^ 
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from which the interface profile was found to be [TP 

V^+V^ = V^. (2) 

This parabolic shape Q describes the non-trivial part of 
the interface where < x^y < t. Outside this region, the 
original boundary is undisturbed. 

Two properties severely constrain the form of possible 
evolution equations for growth inside a three-dimensional 
corner: (a) The governing equation for the interface 
z{x^y;t) must reduce to the two-dimensional form ([T]) on 
the boundaries x = or y = 0; (b) The equation must be 
invariant under the interchange of any coordinate pair. 

Analogously to Eq. ([T]), we seek a three-dimensional 
evolution equation of the form Zt = F{zx^Zy) that in- 
volves only first derivatives (higher-order derivatives are 
asymptotically negligible). The simplest guess is the 
product Zt = [zx/{zx — 1)] [zy/{zy — 1)]. This equation 
reduces to ([T]) on the boundaries x = 0, where Zx = — oo, 
and ^ = 0, where Zy = — oo. The product ansatz is also 
invariant under the exchange x ^ y but not under the 
exchanges x ^ z or y ^ z and therefore is wrong. 

By extensive trial and error, we found that 
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satisfies the necessary coordinate interchange invari- 
ances. These constraints severely limit the form of the 
evolution equation. For example, if we seek a multi- 
plicative correction factor to the product form in (|3| as 
the Laurent series Xl^oo ^n{zx -\- Zy)~^ ^ coordinate inter- 
change invariance gives Aq = l,Ai = —1, while all other 
amplitudes vanish [TQ. Thus Eq. ([3| is the only invari- 
ant choice among the family of evolutionary equations 
parameterized by A^. 

We also found one other elemental evolution equation 
of the form Zt = F{zx-,Zy) that satisfies coordinate in- 
terchange invariance; this form is unique if we again 
seek corrections as a Laurent series representation. This 
second solution is obtained by replacing the factor in 



the square brackets in (|3| with \i ^ {zxZy — Zx — Zy) . 
This equation, which can be re-written more elegantly as 
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and Eq. (|3| are two functionally independent three- 
dimensional evolution equations that satisfy coordinate 
interchange invariance. We believe, but cannot prove, 
that other elemental evolution equations do not exist. 




FIG. 2: (color online) The interface (jsj). 

Our conjecture is that (|3| is the correct evolution equa- 
tion. Evidence in favor of this statement also comes from 
the excellent agreement with simulation data. For this 
comparison, we solve Eq. ([3| by the method of charac- 
teristics. Starting from an empty corner, we find ^ 
that the interface profile admits the following parametric 
representation (Fig. [2]) 
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with q = Zx-, r = Zy and — oo<g, r<0. As a consistency 
check, note that for r = — oo, we have x/t = {q — 1)~^, 
y/t = 0, and z/t = q'^{q — . Eliminating q, we get 
^/x + ^/z = V^, thereby recovering Eq. Q for the inter- 
section of the interface ([5| with the y = plane. 

It seems impossible to eliminate the parameters (g, r) 
from Eq. ([5| and obtain a closed-form representation 
of the interface in terms of x^y^z and t as in the two- 
dimensional case. However, the intersections of the in- 
terface ([5| with certain planes admit simplified descrip- 
tions. For example, for the plane x = y^ corresponding 
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to g = r, we obtain 

X _lz 3 /z\2/3 1 
't ~ 2t~ AVt) ^ 4 ' 

which agrees weh with simulations (Fig. [3| . 
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FIG. 3: (color online) Scaling plot of the surface profile zjt 
versus xjt along the diagonal x = ^ at different times. Upper- 
left inset: the difference A between the simulated values of 
the left- and right-sides of ([6|. Lower-right inset: convergence 
of the diagonal interface speed versus t~°"^^. 



Two additional tests suggest that the conjectured evo- 
lution equation ([3| and its solution ([5| describe corner 
growth accurately. Consider first the advance of the in- 
terface along the ray x = y = z. From ([5|, the position 
of this point is given by [19] 
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Numerically, we measure w ^ 0.1261(2), which agrees 
with our prediction w = 0.125 to within 0.9%. As a sec- 
ond test, we compute the total volume V beneath the 
growing interface at time t. Since the linear dimension of 
the interface grows linearly with time, V = vt^ . To de- 
termine the amplitude we use the parametric solution 
([5| and change from the physical variables (x, y) to the 
parametric coordinates (q^r)^ from which the amplitude 
V reduces to the integral 

J — OO J — oo 



d{q,r) 



We compute the Jacobian ^q^'^^ and the integral using 



Mathematica and find 

3^2 



v = ^ = 0.014457... 



(8) 



Numerically, we measure v ~ 0.01472(3), which is within 
1.8% of our prediction. 

While Eq. ([3| accurately describes the corner inter- 
face, small discrepancies between our measurements of 



the coefficients w and and their predicted values 
(8) persist. The alternative elemental evolution equation 
(4) leads to the interface profile 
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which is the natural generalization of Eq. ([2|. The cor- 
responding values w = ^ and = ^ that arise from 
this profile substantially disagree with simulation results, 
suggesting that Q is wrong. 

From the elemental equations (|3| and Q, we can 
also form two distinct one-parameter families of invariant 
equations p!5|; an additive family 
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and a multiplicative family 
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where for both families the limit c = reduces to (|3| 
and the limit c = — 1 reduces to Q. For the multi- 
plicative class of evolution equations ( |10b| ), the choice 
c ~ 0.074 provides the best fit for the simulated value of 
V [20]. Similarly, for the additive class of equations, the 
optimal mixing parameter is c 0.079. However, a phe- 
nomenon as minimalist as corner interface growth should 
be described by a simple equation that does not contain 
an anomalously small mixing parameter. This aesthetic 
consideration, in conjunction with our numerical results, 
suggest that Eq. (|3| describes corner interface evolution. 

The small discrepancies between our simulation re- 
sults and the predictions that follow from Eq. ([3| (see 
the insets to Fig. |3| suggest that the approach to the 
asymptotic state is slow. A similarly slow convergence to 
asymptotic behavior occurs in various well-understood 
one-dimensional growth models (see e.g. Refs. [21, 22 ). 
For example, for 1+1 dimensional corner growth, the in- 
tersection of the interface with the (1,1) direction evolves 
according to [3| [121 IE] 



(11) 



where ^ is a stationary random variable with (^) > 0. 
Thus averaging over many realizations gives an effective 
velocity w^^^ — \ ^ t~^l'^ . 

For growth inside a three-dimensional corner, we there- 
fore anticipate that w^^ ~ \ ^ ^ with a still-unknown 
exponent a. Very extensive simulations for flat interfaces 
in 2+1 dimensions indicate that a is close to 0.77 [23H25] . 
On the other hand, extrapolation from our simulations 
for t < 20000 suggests that a ^ 0.74. This difference 
in exponent estimates suggests that t = 20000 is still 
outside the long-time regime for growth inside a three- 
dimensional corner. This slow approach to asymptotic 
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behavior could be the source of the discrepancy between 
our simulation results and the theoretical prediction (|3| 
for the interface profile. 

Our argument for the form of the evolution equation 
can be generalized to higher dimensions. Applying coor- 
dinate interchange invariance and related symmetry con- 
siderations, we conjecture that in d dimensions the height 
h{xi^ . . . , Xd-i]t) satisfies 
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where hi = ^ . These equations are again solvable using 
the method of characteristics [15 . 

We emphasize that computing the limiting shape — 
the primary characteristic of the interface — represents 
only a first step to understanding its properties. One 
challenging problem, given that interface fluctuations are 



unknown even for flat interfaces, is to generalize Eq. (11 ) 



to account for fluctuations of an interface that grows at 
a three-dimensional corner. Also of interest are height- 
height correlations at different locations and different 
times. In 1+1 dimensions, these correlations decay slowly 
along the characteristic curves of the evolution equation 
[7l[26]. Whether similar behavior occurs in 2+1 dimen- 
sional corner growth is unknown. 

Fluctuations of integral characteristics of the interface, 
such as the crystal volume, may be more tractable and 
give rise to new phenomena. Consider, for example, the 
total number of sites of various fixed degrees (number of 
adjacent vertices). Sites of degree 3, in particular, can be 
categorized as either inner or outer corners. The num- 
ber of inner corners grows as A^in = ^ = 3i't^, with 
V = 37r^/2^^ to leading order. One might anticipate the 
same asymptotic growth for outer corners, but simula- 
tions indicate that the latter grows slightly faster [15 : 
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Note that in two dimensions Ain — Aout = 1- For Ising 
corner growth in three dimensions, Ni^ — A'out is also 
positive and grows with time as t^/^. This makes the 



behavior in ( 13 ) quite puzzling. 



The other major challenges are to generalize from strict 
corner growth to Ising growth, where adsorption at in- 
ner corners and desorption from outer corners occur with 
equal rates, and to equilibrium interfaces, where the des- 
orption rate exceeds the adsorption rate. The corre- 
sponding equilibrium shape has been determined both 
in two [27 and three dimensions [28 , and its shape fluc- 
tuations have also been studied [29] . In analogy with the 



conjectured evolution equations (12) for corner growth. 



there may also exist an exact generalization of equilib- 
rium limiting shapes [271 EH] in higher dimensions. 
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